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1. Introduction 

Quantum graphs provide models in many disparate areas of mathematical physics. They 
have been used to describe electrons in organic molecules, nanotechnology, wave guides, 
Anderson localization, quantum chaos and mesoscopic systems to name but a few, see 
[35] for a review. One of the reasons why quantum graphs have been so successful is 
that they possess an elegant exact trace formula expansion for the spectrum in terms 
of periodic orbits on the graph. Such formulae were first derived by Roth [39] and 
independently by Kottos and Smilansky [331 EI], see [11] for the current state of the 
art. The prototype for quantum mechanical trace formulae was provided by Gutzwiller 
who showed that in the semiclassical limit the energy spectrum can be expressed as a 
sum involving an infinite number of classical periodic orbits [19J. For most systems this 
trace formula suffers from convergence issues that can only be resolved by limiting 
the contribution of long periodic orbits and in the process valuable information is 
suppressed. One way to recover this information is through a remarkable resummation 
procedure which connects long and short pseudo orbits (collections of periodic orbits) 
via a truncated quantum zeta function. The resulting expression, discovered by Berry 
and Keating [HI [TJ, [28] and derived using a different method by Bogomolny [9], is termed 
the Riemann-Siegel lookalike formula due to the connection with its number theoretical 
namesake. It allows one to approximately determine the quantum spectrum up to some 
given energy, using only a finite number of periodic orbits. 

For quantum graphs, there exists both an exact Riemann-Siegel lookalike formula 
and an exact pseudo orbit expansion including only a finite number of pseudo orbits 
with at most as many bonds as the graph. Kottos and Smilansky describe this without 
an explicit formula in [M]. Here we provide a basic argument to obtain a minimal 
pseudo orbit expansion. This also sheds light on the origin of the cancelations seen 
in such pseudo orbit expansions. We use the expansion to examine the coefficient of 
the characteristic polynomial and to produce a simple expression for a secular equation 
whose roots provide the spectrum of the graph ^ Ai ^ A2 ^ . . .. This enables the 
whole spectrum of the graph to be formulated in terms of a finite number of pseudo 
orbits which visit at most the same number of bonds as are in the graph. 

To demonstrate the efficacy of this minimal pseudo orbit expansion we describe two 
applications to the spectral determinant and vacuum energy of graphs. Formally the 
spectral determinant is the product of the eigenvalues, 



where A is some spectral parameter and the prime indicates that zero modes are omitted 
where present. The spectral determinant on graphs has been investigated by a number 
of authors [31 H21 [131 HH HS1 [231 EH [261 [29] ■ In particular in pQ Akkermans et al. 
derive a formula for the spectral determinant of the graph as a product over all periodic 
orbits. After expanding the product they discuss cancelations between certain sets of 
periodic orbits using a diagrammatic approach. While the pseudo orbit expansion that 




(1) 



n=l 
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we describe in section IfTTI has not been written down previously, it can be obtained using 
the cancelation mechanisms from p] , as is explained in the appendix. The method that 
we employ in the body of the paper is more straightforward and the proof shows that 
the existence of such a pseudo orbit expansion does not depend on the structure of the 
scattering matrix, as implied in [TJ. 

A second application of the pseudo orbit expansion allows us to provide a trace 
formula over the short pseudo orbits for the vacuum energy of the graph. Formally the 
vacuum energy is 

^jE'v^- ( 2 ) 

71=1 

Vacuum energy is responsible for the Casimir effect, conventionally seen as an attraction 
between uncharged parallel conducting plates at short distances. Analytic results for the 
vacuum energy in most geometries are hard to obtain and consequently it is interesting to 
consider mathematical models of vacuum energy where the analysis can be taken further 
[38] . In [17] Fulling et. al. consider the vacuum energy of a star graph with Neumann 
conditions at the vertices and show that the Casimir effect is repulsive. The vacuum 
energy for generic graphs and vertex conditions was investigated in [2j [23J EU [25] . 

To regularize both the spectral determinant and vacuum energy we employ the 
spectral zeta function, 

oo 

C(s,A) = ^'(A + A n r, (3) 

71=1 

which is convergent for 3?s > 1. Making an analytic continuation to the complex s plane 
the spectral determinant is defined to be, 

S(A):=exp(-C'(0,A)). (4) 

Similarly the regularized vacuum energy is obtain from the £(—1/2, 0), see section [6~2l 
The article is laid out as follows. Section [2] introduces quantum graphs and the 
scattering approach used to define a secular function. We then provide, in section El a 
simple derivation for the pseudo orbit expansion of the secular function. The following 
sections are devoted to applications of this expansion. Section H] examines the variance of 
the characteristic polynomial coefficients, section [5] develops a regularized expression for 
the zeta function and section [6] uses this expression to obtain pseudo orbit expansions for 
the spectral determinant and vacuum energy. After summarizing the results we include 
an appendix which describes alternative cancelation mechanisms which provide pseudo 
orbits expansions. 

2. Quantum graphs 

Let Q = (V, B) be a graph with the finite sets of vertices V and bonds (or edges) B. 
Each bond, b G B, connects a pair of vertices, u, v e V, see figure [TJ We use the notation 
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b = (u, v) to denote a directed bond where o (b) — u and t (b) = v, with o and t standing 
for the origin and the terminus of b. 

We denote by b = (v, u) the reverse of the directed bond b. For every bond beB 
we insist that b is also in B. The size of the set B is then denoted by IB := \B\. The set 
of all bonds whose origin is the vertex v is denoted B v and its size is called the degree 
of the vertex, d v := \B V \. 

The combinatorial graph Q becomes a metric graph G when we identify the bonds 
of Q with one dimensional intervals [0, Z&] with a positive finite length lb {l b = lb). We 
can then assign to each bond b a coordinate, Xb, which measures the distance along the 
bond starting from o(b). The reverse directed bond, b, posses the coordinate, x b , and 
the two coordinates are related by Xb + x b = l b . We denote a coordinate by x, when its 
precise nature is unimportant. 2£ = X^ee L b is twice the total length of G. 




Figure 1. (a) A fully connected graph with four vertices and six bonds and (b) a graph 
with five vertices and eight bonds. 



The quantum graphs we consider are metric graphs equipped with a self-adjoint 
differential operator ri, the Hamiltonian. Here we are particularly interested in the 
negative Laplace operator, 

ri : /(*) ^ -0 , (5) 
or the more general Schrodinger operator, 

U : fix) H- -0 + V(x)f(x) , (6) 

where V(x) is a potential, which we assume to be bounded and piecewise continuous. 
Note that the value of a function or the second derivative of a function at a point on the 
bond is well-defined, thus it is not important which coordinate, Xb or x b is used. This 
is in contrast to the first derivative which changes sign according to the direction of the 
chosen coordinate. 

To complete the definition of the operator we need to specify its domain. We denote 
by H 2 (G) the space 

i7 2 (G):=0ff 2 ([O,4]), (7) 

6=1 
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which consists of the functions /onG that on each bond b belongs to the Sobolev space 
H 2 ([0, l b ]). The restriction of / to the bond b is denoted by where fb(xb) = fi(h~ x b)- 
The norm in the space H 2 (G) is 

B 

\\f\\H*{G) ■= 1 1 ^11 Wad • ( 8 ) 

6=1 

Note that in the definition of H 2 (G) the smoothness is enforced along bonds only, 
without any matching conditions at the vertices at all. However, the standard Sobolev 
trace theorem (see e.g. [15]) implies that each function fj, and its first derivative have 
well-defined values at the endpoints of the bond b. 

All conditions matching functions at the vertices that lead to the operator ([5]) 
being self-adjoint have been classified in j22j E21 [36] . It can be shown that under these 
conditions the operator "H is bounded from below |36j . In addition, since we only consider 
compact graphs, the spectrum is real, discrete and with no accumulation points. Thus 
we can number the eigenvalues in the ascending order and denote them with {\ n }™ =1 . 
The n'th eigenvalue obeys the following equation 

d 2 f n 



dx 



+ V(x)f n (x) = X n f n (x) , (9) 



where f n is the corresponding eigenf unction. We also use k n , such that A n = k 2 , and 
say that {k n }™ =1 is the /c-spectrum of the graph. 

In this paper we will primarily consider graphs where the domain of the operator 
([6]) consists of the functions / G H 2 (G) which obey Neumann conditions at all vertices 
of the graph. The function / is continuous at the vertex v, i.e. /^(O) = /b 2 (0) for all 
bonds bi, b 2 £ B v . The derivatives of / at the vertex v satisfy, 

E^(o) = o. do) 

We will assume that G contains no loops (bonds b with o(b) = t(b)) and no pairs of 
vertices connected by multiple bonds. However, as a loop or bond can be subdivided by 
introducing a vertex of degree two with Neumann boundary conditions without changing 
the spectrum, this assumption does not reduce the generality of the results. 

In addition, we will assume that the graph's spectrum is non-negative, Ai > 0. This 
assumption holds in the case of graphs with Neumann conditions with no potential, 
where Ai = 0. The condition can be relaxed by accepting some additional technicalities, 
see [31]. The next section describes an approach to obtain the spectrum of such a graph. 

2.1. The scattering approach to the graph spectrum 

Let G be a graph with Neumann vertex conditions and no potential. The eigenvalue 
equation, 
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on each bond has a solution that can be written as a linear combination of two complex 
exponentials if k ^ 0. We write 

f b (x b ) = a b a e- ikx " + a° b ut e [kx » , (12) 

for the solution on the bond b. The same function can be expressed using the coordinate 

f b (x h ) = af e- ik *> + af ft . (13) 

Since these two expressions should define the same function and since the two 
coordinates are related, through the identity x b + x b = l b , we arrive at the following 
relations 

fl in = e iW 6a out &nd fl m = e ikl ba out _ ^ 

Fixing a vertex v and using the matching conditions to relate solutions /b for all 
bonds b with o (b) — v one arrives at 

sz* = *to(k)a*, (is) 

where a° ut and a™ are vectors of the outgoing and incoming coefficients at v and a^ v \k) 
is a d v x d v unitary matrix, d v being the degree of the vertex v. The matrix cr^(fc) 
is called the vertex-scattering matrix. It is k independent for the Neumann vertex 
conditions and its entries were calculated in [34J: 

= I" - V • (16) 

Collecting all coefficients a™ on the whole graph into a vector a of size IB such 
that the first B entries correspond to bonds which are the inverses of the last B entries. 
We can then define the matrix J acting on a by requiring that it exchanges a™ and a™ 
for all b such that, 



J= • (17) 




Then, collecting equations (115]) for all vertices into one system and using ffl4|) we have, 

Je~ ikL a = E{k)a , (18) 

where L = diag{7i, . . . , Ib, h, ■ ■ ■ , Ib} is a diagonal matrix of bond lengths and S (fc) is 
block-diagonalizable with individual a^ v ' (k) as blocks. This can be rewritten as (note 
that J- 1 = J), 

a = e ikL JZ{k)a . (19) 

The unitary matrix S(k) := JS(/c) is called the bond-scattering matrix. The unitary 
matrix U (k) := e lkL S(k) is the quantum evolution operator. Using this notation we 
write (fT9~j) as 

(I- U(k))a = . (20) 
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A vector a satisfying equation (120 j) defines an eigenfunction on the graph with the 
eigenvalue E = k 2 . Hence we arrive at the conclusion that the non zero eigenvalues of 
the graph are the solutions of 

det (I—U (k)) — . (21) 

This is called the secular equation. It is the major tool for investigating various 
spectral properties of quantum graphs. The interested reader can find a more elaborate 
description of the derivation above in [18]. We end by noting that a similar equation 
may be derived for the case of a potential which differs from zero [3D] . 



3. A pseudo orbits expansion for the secular function 

We define the secular function 

F (k) := det (U (k))~^ det (I - U (k)) . (22) 

This is the LHS of the secular equation ( !2~TI) . multiplied by det (U (k))~^ to turn it into 
a real function for all k G IR, as U (k) is unitary. We recall that all positive values of 
the graph's fc-spectrum are given by the zeros of F (k). 
The characteristic polynomial of U (k) is, 

2B 

F € (k) := det (£1 - U (A;)) = E a n e B ~ n , (23) 

n=0 

where the fc- dependence of the coefficients is omitted for brevity. One can express the 
secular function in terms of the coefficients of the characteristic polynomial 

F (k) = det (U F 1 (k) 

= a 2B ° n ' ( 24 ) 

ra=0 

The main result of this section is a theorem which expands the secular function as a 
sum over pseudo orbits on the graph. Before stating the theorem, we describe orbits on 
graphs. 

A periodic orbit on G is a closed path that starts and ends at the same vertex. 
We denote such a trajectory by the corresponding ordered set of its bonds, 7 := 
(61, &2> • • • b n ), such that Vi, t(bi) = o{b i+ i) and t (b n ) = o(b\). Cyclic permutations 
of the bonds of 7 are the same periodic orbit, for example, 7 = (b 2 , . . . b n , 61). 

A primitive periodic orbit is a periodic orbit which cannot be written as a repetition 
of a shorter periodic orbit. We can thus define the repetition number of an orbit r 7 as 
the number of repetitions of the primitive periodic orbit which it contains. For example, 
if 70 = (bi, b 2 , 63) is a primitive periodic orbit, then 7 X = (pi, b 2 , b 3 , bi, b 2 , 63) is a periodic 
orbit with repetition number r 71 =2. Note that the repetition number of a primitive 
periodic orbit is always one. For a periodic orbit 7 we use 5 7 := n to denote the number 
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of bonds in 7, the topological length of 7. Similarly the metric length (or length) of 7 is 
denoted Z 7 := XVe 7 ^he P r °duct of scattering amplitudes around 7 is; 

v4 7 = Sb 2 b 1 Sb a b 2 ■ ■ ■ Sb„b„- 1 Sb 1 b n , (25) 

where we have suppressed the k dependence if present. 

We now define a more complex structure which we call a pseudo orbit (also called 
a composite orbit in the literature). A pseudo orbit is a collection of periodic orbits, 
7 = {71,72, • • -7m}- An irreducible pseudo orbit 7 is a pseudo orbit which contains no 
directed bond more than once. So every periodic orbit jj in 7 has at most one copy of 
each directed bond and every pair, 7^, 7^, has no directed bond in common. In particular 
this implies that an irreducible pseudo orbit consists only of primitive periodic orbits 
and the primitive periodic orbits are distinct. 

We will use the following notation for properties of 7: 

the number of periodic orbits in 7. (26) 
the total topological length of 7. (27) 

the total metric length of 7. (28) 

(29) 

The null pseudo orbit is the pseudo orbit which contains no orbits. We treat it as an 
irreducible pseudo orbit with the following properties: mg = Bq = Iq = and Aq — 1. 

Theorem 1. The coefficients of the characteristic polynomial F^ (k) are given by 

a n = (-1)"^ ^7 (AO exp (ikl^) . (30) 

Proof. We first rewrite the characteristic polynomial using the permutation expansion 
of the determinant 

2B 

det (£1 - U (*)) = s S n (P) U^ I U (*) > ( 31 ) 

p£S 2B 6=1 

where the indices of the matrix, p (b) , b, label directed bonds of Q. We may now use the 
unique presentation of the permutations as a product of disjoint cycles 

P = (h, h,--- Ki) (K1+1, ■ ■ ■ , b ni +n 2 ) ■ ■ ■ y } J2 m =C 1 n i+ v ' ' ' ' ^ 

for a permutation p with m p cycles where the j'th cycle has length rij. The notation 
above means p acts on the 2B directed bonds such that p(&i) = &2 etc. The 
representation (1321) is unique up to cyclic permutations of each cycle and no directed 
bond appears more than once. 



m 7 :- 
.By := 



M 

E' 

7j£7 

n-4 



7j67 
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Recall that [U (k)] , h \ b 7^ only if the directed bond p (b) is connected to the 
directed bond 6, i.e., o(p[b)) = £(6)jj] From which we conclude that a permutation 
p ^ I can have a non zero contribution to the sum ( 13TI) only if o(p(b)) = t (b) for all 
b in p. This allows us to interpret (|32|) as an irreducible pseudo orbit on the graph. 
Each periodic orbit in the pseudo orbit corresponds to a cycle of p. The pseudo orbit 
is irreducible since no directed bond appears more than once in (1321) . The fixed points 
of p, the bonds b not in p so b = p (b) , each produce a multiplicative factor of £ in the 
contribution of p in ( 13TT) . 

The coefficient a n multiplies £ 2B_n and is therefore obtained as a sum of all 
permutations p G S2B with 21? — n fixed points. Each such permutation corresponds 
to an irreducible pseudo orbit, 7, whose bond length is = n. Using the notation of 
( 1321) . the pseudo orbit p consists of m p orbits and n = Y^7=x n i- The contribution of 
this permutation (irreducible pseudo orbit) to ( 13~TT) is given by 

sgn(pK 2S -"(-l)^^ (U b2bl U b3h ...U blbni ) (U bni+2bni+1 ...U bni+lbni+n2 )..., 

where we omitted the k dependence from U for clarity. If we recall that U b2bl = 
exp (ikl b2 ) iSfegbj., we see that the contribution above (up to a sign) is £ 2B_n exp (ikly) Ay. 
The result ( 130]) now follows since 

sgn(p) (-1)^7=1^ = (_i)Er="iK+i) (_i)Er= P !%' = (_i)"p = (_!)-p . 

□ 

We now employ the unitarity of U (k) to get a Riemann-Siegel lookalike formula 
which connects pairs of coefficients of the characteristic polynomial [34J. The 
characteristic polynomial can be rewritten 

F € (k) = det {-iU (k)) det (^^I - U ] (k)) 

2B 

= det(-ZU(k))Y,<C- 2B 

n=0 

2B 

= det (U(k))J2<C (33) 

n=0 

Comparing this to (k) = Y^^=o a ni 2B ~ n an d using a 2 s = det (U (k)), we deduce the 
following symmetry relation of the characteristic polynomial coefficients 

a n = a 2B a?2B-n ■ ( 34 ) 

This symmetric expression allows us to express the secular function in terms of shorter 
pseudo orbits, whose total number of bonds is less than or equal to the number of bonds 
of the graph. 

f Provided that all elements of the vertex scattering matrices are non-zero. This would not hold, for 
example, in the case of equi-transmitting matrices |27j where back-scattering is prohibited. 



Finite pseudo orbit expansions for spectral quantities of quantum graphs 10 

Theorem 2. The secular function has the following expansion in terms of irreducible 
pseudo orbits, 

F{k) = 2 (-1)™" ^ 7 cos(£;(Z 7 -£) - 6 (k)) H (B - By) , (35) 

l\B^<B 

where det S = exp {¥26 (k)), and H is the Heaviside function, see $39\). 



Proof. Plugging the symmetry relation f[34|) in the secular function expression (j24j) we 
obtain, 

25 

F{k) = det(U(k))~* J2 an 



n=0 

5-1 



— ^2 ( fl 2^ fln ~*~ fl 2B a n) + a 2B a B 
n=0 

f ^ 1 

= 3? < 2 a 2jB 2 a n + a 2 J a B ? , (36) 

I n=0 J 

— 1/2 

where a 2S a# is a real number which can be seen from the symmetry relation as 
1 025 1 = 1- Combining this with theorem [1] we see 

f 

F(k) = ft I 2(detS)-^J2 J2 (-1)™" ^exp(ikm - £)) 

^ n=0 7I B^=n 

+ (detS)~5 (-l) m "^7exp(iA;(^-£)) I . (37) 

7|5,=5 J 

The unitarity of S allows to denote det S = exp (i29 (k)) and to rewrite the above as 

6-1 

F(A;)=2^ ^ (-l) m - A 7 cos (A; (Z 7 -£)-#(£;)) 

n=0 7 | B^=n 

+ ^ (-l) m ^ 7 cos(A;(/ 7 -£) -0(Jfe)). (38) 

7|5,=B 

This expression is made more compact by using the Heaviside step function 

{0 x < 
\ x = . (39) 
1 x > 

□ 

Finally, if we make the assumption that our graph has Neumann conditions at the 
vertices then either 9 = or 9 = | with no /c-dependence. This results in an expansion 
in terms of either cosine or sine functions. In the remainder of the article we will assume 
that the graph has Neumann vertex conditions. 
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We also note that the point of truncation < B in (1351) is analogous to a cutoff 
at half of the Heisenberg time. The Heisenberg time T# is 2tt multiplied by the mean 
spectral density, which for a graph is C/tt, hence Th = 2C If Z av is the average length 
of a bond on our graph, and we assume that the bond lengths are all close to / av , then 
ps l av B^ and Th/2 = l av B. Thus the cutoff B^ < B resembles < Th/2 used in the 
semiclassical expansions [SI U\ [28] . This in turn originated with the cutoff introduced 
in the functional equation of the Riemann zeta function, known as the Riemann-Siegel 
formula. However, the argument leading to the cutoff in ( |35i) is due to the unitarity of 
U(k). A similar approach applies to the dynamical zeta function [HI [9]. 



4. Variance of coefficients of the characteristic polynomial 

As a first application of these results we are now in a position to write a simple 
expansion for the variance of coefficients of the characteristic polynomial ( l30i) . The mean 
values of the coefficients (other than = 1) are zero as the average over k of e lkl ^ is 
necessarily zero for pseudo orbits of topological length n > 1. The second moment of the 
coefficients of the characteristic polynomial, the variance, was investigated numerically 
in [31] where it is seen that they may deviate from the random matrix predictions 
even for graphs whose spectral-statistics match those of appropriate random matrix 
distributions. Tanner [12] considered the variance of the coefficients for a related graph 
model which generates a unitary stochastic matrix ensemble [H] over which the average 
is taken. 

From theorem [1] the variance of the coefficients of the characteristic polynomial can 
be expressed as a sum over pairs of irreducible pseudo orbits 7,7', 

1 r K 

{\a n \ 2 )k = (-l) m ^ +m ^'A^ lim — / e ik ^- l ^dk 



= (-1)^^^,!,' (40) 

So a pair of irreducible pseudo orbits of topological length n only contributes to the 
sum when the total metric lengths of the pseudo orbits are equal. When the set of 
bond lengths is incommensurate this implies that the two irreducible pseudo orbits visit 
a given pair of directed bonds {b, b} the same number of times. Hence the variance 
is independent of the bond lengths, provided that they are incommensurate. The 
formulation of the variance ( f4"0l) also holds for vertex scattering matrices that do not 
correspond to Neumann matching conditions provided that the scattering matrices are 
independent of k. 



4-1. Diagonal approximation 

In the case of incommensurate bond lengths the first order contribution to the variance 
will come from pairing an irreducible pseudo orbit 7 with itself, 7' = J, or with any 
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other orbit 7' obtained from 7 by reversing the direction of a subset of orbits in the 
pseudo orbit, as in the diagonal approximation of Berry [5]. As 7 contains distinct 
orbits there are typically 2 m<? different pseudo orbits that can be obtained from 7 by 
reversing a subset of the orbits in 7- Making this diagonal approximation we obtain 

(\a n \ 2 ) k ^ 2 ^4- ( 41 ) 

j\Bj=n 

Note that while making the diagonal approximation, we did not take into account the 
existence of self retracing orbits. These orbits are equal to their reversed orbit. Pseudo 
orbits which contain self retracing orbits should therefore have a smaller prefactor in 
f )4T|) . However, the proportion of self retracing orbits among all periodic orbits decreases 
rapidly with their length, which justifies the approximation. 

As A? < 1 this suggests that the dominant contribution in the diagonal 
approximation will come from pseudo orbits comprised of many distinct periodic orbits 
where is maximal. Figures |2] and [3] show the variance of the coefficients of the 
characteristic polynomial for the two graphs shown in figure CD computed from the 
pseudo orbit expansion (HOjl , the diagonal approximation ( 14T1) and by directly averaging 
coefficients of (k) over a range of approximately 2000 mean level spacings. The 
diagonal approximation works well for these examples. One can also see that where 
the diagonal approximation fails, around n/2B = 1/2, the diagonal approximation over 
estimates the actual value which suggests that for these small graphs the over counting 
of the self-retracing orbits is significant. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



n/2B 

Figure 2. Variance of coefficients of the characteristic polynomial of a fully connected 
graph with four vertices (shown in figure W(a)) with random bond lengths chosen 
uniformly in [0,1]. Numerical results were averaged over 200 000 values of k over 
a range of approximately 2000 mean level spacings and are compared with the pseudo 
orbit expansion, diagonal approximation and random matrix theory. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



n/2B 

Figure 3. Variance of coefficients of the characteristic polynomial for the graph shown 
in figure{l](b) with random bond lengths chosen uniformly in [0, 1]. Numerical results 
were averaged over 200 000 values of k over a range of approximately 1700 mean level 
spacings and are compared with the pseudo orbit expansion, diagonal approximation 
and random matrix theory. 



The variance of the coefficients of the characteristic polynomial of a 2B x IB 
random matrix were obtained in [21]. For random matrices from the Circular Orthogonal 
Ensemble (COE) and Circular Unitary Ensemble (CUE) (/3 = 1 and (3 = 2 respectively) 
the variance of the coefficients are given by, 

f 1 + n{2B~n) R _ , 

(K\ 2 ),= l [ + 2B+1 P p J 2 ■ (42) 

In the semiclassical limit the spectrum of generic quantum graphs is seen to 
approach that of an appropriate random matrix ensembles. This goes back to the 
results of Kottos and Smilansky for quantum graphs [31] which is the graph setting of 
the well known correspondence of Bohigas, Giannoni and Schmit [10j. The semiclassical 
limit is the limit of increasing spectral density, which for quantum graphs is the limit of 
a sequence of graphs with increasing number of bonds [18]. In [I] Berkolaiko and Winn 
formalize the connection between the spectrum of a scattering matrix and the spectrum 
of the corresponding Laplace operator on the graph. Of particular relevance to the 
current discussion, when the spectrum of the graph is distributed like the eigenvalues of 
a random matrix we also expect that the eigenphases of the scattering matrix on the unit 
circle will have the same distribution. When this is the case one must also anticipate that 
the coefficients of the characteristic polynomial of the graph are distributed according to 
the random matrix prediction. As this was not observed it appears that the variance of 
the coefficients of the characteristic polynomial is more sensitive to the size of the graph 
used for the numerics than the spectrum itself. In [12] Tanner investigated sequences 
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of binary graphs without time-reversal symmetry where he observes converges to a 
constant different from the random matrix result and dependent on the classical limit 
of an underlying Markov process. 



5. Zeta function 



Let z = k + it with k,t real, following theorem [2] we will denote 

F{z) = 2 (-l) m ^7 cos (^7- C)-6)H(B- By) (43) 

*/\B^<B 

so that the secular equation reads F(z) = 0. F(k) is an even or odd function of k 
depending on whether 9 is or 7r/2, so we have a sum of either cosine or sine functions 
respectively. Notice that F(k) must be even or odd as the fc-spectrum is symmetric 
about zero. For Neumann vertex conditions in the limit z — > 0, 



F(z) 



~ z 

2->0 



R-V+2 



+ 0{z 



B-V+4\ 



which comes from the identity dim{ker(I — {7(0))} = B — V + 2 [37J. We define 



z 



B-V+2 



(44) 



(45) 



Importantly F also has roots k n but does not vanish at zero, F(z) ~ 1 + 0(z 2 ) as 
z — > 0. We see that F(—z) = F{z). To formulate the spectral zeta function we follow a 
procedure introduced for graphs in [21] which we only sketch here. We take the function 
F(z) and use the argument principle, evaluating an integral around a contour c enclosing 
the positive real axis, see figure BJa). 

d 



C(«,A) = -ir l(z 2 + \y 

2m 

which converges for 3?s > 1. 

(a) .. 
it 



dz 



log F(z) dz , 



(46) 



-iVA 



(b) 



it 



( X x - 



-x — x- 



x zero of F(z) 



— x^ X- 



zero of F(z) 



k 



Figure 4. Contours used to formulate £(s, A) (a^) before, and (b) after, the contour 
transformation. Branch cuts are represented by hashed lines. 
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C(s, A) = / (t 2 - X)- s - log F(it) dt (47) 

7T 



We may deform the contour c to d where the integral is along the imaginary axis 
figure 0Jb). 

/•OO 1 

sm7rs / ,. n . d 
JVx dt 

The segment between ivA and — i^/A does not contribute to the integral due to the 
symmetry of F(z). Analysis of the integral shows that the representation of the zeta 
function ( l47j) converges in the strip 1/2 < 3?s < 1. As we will see, the restriction 
> 1/2 is due to the t — > oo behavior of F(it). 

= E (-!) m ^7 ± e* (/ ^ £) ) - B,) , (48) 

where the choice of sign for e 1 ^ - ^ depends on whether = or 7r/2q] As no directed 
bond appears in 7 more than once < C. Hence 

which is the contribution of the null pseudo orbit 0; the pseudo orbit of length zero. 
Consequently, 

^logF(it) ~ C~ B ~ V + 2 , (50) 

and we see that (j4"Tj) converges for 9fts > 1/2. To make the analytic continuation of 
the zeta function to the left of the line 9fts = 1/2 we add and subtract the t — > 00 
asymptotics in the integral. We thus obtain an expression for ((s, A) valid for 3Rs < 1, 

C(s, A) = ^ f> - log(F(it)t™e-«0 dt 

+ C T ±fl§^-- B -^±lx-. (5!) 

2^7TL [SJ 2 

6. Spectral determinants and vacuum energy 

We are now in a position to use the zeta function to formulate the spectral determinant 
and vacuum energy of the graph as sums over a finite set of pseudo orbits. 

6.1. Spectral determinant 

The spectral determinant, formally the product of the shifted eigenvalues (JTJ), was 
defined as <S(A) = exp(— £'(0, A)), see [24J. From ( 15T|) we see 

C'(0, A) = [log(F(it)t B ^ +2 e-^)] ~ - CVX + B ~ V 2 +2 log A 

= log - log(#(h/A) A^e-^) - CVX + - - j + 2 log A . (52) 

f The sign of the whole function is not important as multiplying F by a constant factor does not 
change its properties. 
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Hence the spectral determinant is given as 
2 

xi B-V+2)/2 



5 ( A ) = x(B-L)/2 E (-!) m ^7 cosh(VA(^ - £)) B(B - By) , (53) 



when = or the same formula with cosh replaced by sinh when 9 = it/ 2. This 
concisely formulates the spectral determinant of the graph in terms of a finite set of 
dynamical quantities. 

6.2. Vacuum energy 

A second application of the pseudo orbit approach is the zeta-regularized vacuum energy. 
The vacuum energy, formally the sum of the square roots of the eigenvalues §Z§, is defined 

as 

£ C =±C(-V2,0) . (54) 

Were C(s, 0) to have pole at s = —1/2 a more general definition would be required, see 
e.g. [30]; however, we will see that there is no pole in the present case. The Casimir 



force on a bond b is given by — -^£c- 



The zeta function for A = is 



sin tts f°° ._ 2s d 



C(s, 0) = / r 2s - log F(it) dt , (55) 

7T J dt 



in the strip 1/2 < 3?s < 1. Splitting the integral at t — 1 

C(S) 0) = ^^\- 2s ^logF(it)dt + ^r 2 ^i ogj F(it)d^ , (56) 

and the restriction to 5Rs > 1/2 comes from the t — > oo behavior in the second integral. 
Adding and subtracting the asymptotic behavior (I49p in the second integral (as in (1511) ) 
we obtain, 

sin7rs f 1 9c d , - .. s , sin7rs f°° 9o d , , R T/ , 9 , r . , 

C(s, 0) = / r 2s - log F it dt + / t- 2s - \og(F(it)t B - v+2 e- tc ) dt 

7T J dt 7T J 1 dt 

sinvrs / £ V - B - 2\ 

+ hr— r + ^ > (57) 



7T V 25 - 1 2S 

which is now valid for 3?s < 1. Using this representation we obtain 
£ c = / t— log F(it) dt + — + 



2?r ,/ dt 4?r 2tt 

1 v-^ f°° tl-e~ tls < - t(2C - l-)e~ t(2C ~ 1 ^ 

+ ^_2^J- 1 ^J 1 E^< B (-l)^A(e^+e-^)) dt - (58) 

This is an exact trace formula for the vacuum energy where the lengths of orbits included 
in the sum are truncated at a finite value. For a particular graph the integrals can be 
evaluated numerically. 
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7. Conclusions 

We have investigated a pseudo orbit approach to spectral questions on quantum graphs. 
The possibility of a pseudo orbit expansion of the secular function of a graph was 
described in [34]. Here, however, we have presented the explicit formula. As quantum 
graphs provide a system where the spectral analysis can be carried out exactly, they 
serve as an ideal example in which to investigate the various approaches to cancelations 
which appear in pseudo orbit expansions. To this end while the permutation expansion 
of the determinant provides a neat direct route to the expansion using irreducible pseudo 
orbits, we also describe other methods of deriving this formula in the appendix. In these 
cases a sum over pseudo orbits is seen to contain substantial cancelations when pseudo 
orbits are grouped appropriately. 

Having obtained the expansion of the secular function in terms of pseudo orbits we 
demonstrated how this can be applied to three spectral quantities, the variance of the 
coefficients of the characteristic polynomial, the spectral determinant and the vacuum 
energy. We provide an exact pseudo orbit expansion for these quantities using only 
dynamical information for a finite set of short periodic orbits whose topological length 
is at most the length of the graph. In the case of the vacuum energy it is particularly 
interesting to note that such a quantity can be obtained from a finite set of dynamical 
properties. In more complex physically realistic systems it is often hard to calculate 
the vacuum energy. If it is possible to extend these techniques to compute the vacuum 
energy of higher dimensional systems one may be able to describe the vacuum energy 
precisely using orbits with a finite cut off by length. This would seem to offer the 
possibility of precise numerical evaluation of the vacuum energy of general systems. It 
is a hard problem to provide an accurate numerical evaluation of the vacuum energy, 
such calculations are often based on the proximity force approximation where a rigorous 
justification is lacking. 
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Appendix - More on cancelation mechanisms 

The appendix is dedicated to a description of the various canceling mechanisms which 
play a role in the derivation of pseudo orbits expansions. A combination of some of these 
schemes provides an alternative proof of theorem [TJ Some cancelation mechanisms were 
already presented in the literature in various contexts (not necessarily for quantum 
graphs), as we will mention. However, not all of them were fully developed to present 
the current pseudo orbits expansion. The purpose of bringing them together here is to 
give a state of the art description of the various mechanisms and their interrelations. 




Appendix A.l. From an infinite expansion to a finite one 

We start from the characteristic polynomial of U (k) and use the matrix identity 
det (I - U{k)) = exp Tr In (I - U{k)) = exp Tr (- ^f- 

too 

where we use the notation Tj := TrfC/^A;)^]. 

The traces Tj can be expressed in terms of periodic orbits 

T 3 = fV 7 *, ( A -2) 

llB^j 7 

see notation in section [3j Inserting f ]A.2j) in f lA.lj) we see, 

-EE r i ^ L ) • (a.3) 

i=l 7l B 1= j 7 / 

Fixing a primitive periodic orbit tt and summing over its repetitions, 

/ _ °° Ar irkl^ 

F,(k) = e B ew - Yl £r^^- 



7r primitive r=l 



i 2B exp ( ln (! - C B *Aj kl ") 

\7r primitive , 

e B n t 1 - r^Ar^) . 



(A.4) 



7r primitive 

Expanding the last line above, one would get an infinite sum, each of whose terms is the 
contribution of some collection of primitive periodic orbits. Namely, each term would 
correspond to a pseudo orbit which consists of only primitive orbits; a primitive pseudo 
orbit. This is not to be confused with an irreducible pseudo orbit see section [3] - which 
form a subset of the primitive ones. Keeping in mind that the characteristic polynomial 
has only finite powers of £, F^ (k) = £ 2B Y^Lo a n£,~ r \ we see that the contribution of 

f This holds for k + ie which is to be understood in the limit e —> 0. 
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many of these pseudo orbits cancel. The remaining pseudo orbits give the following 
expression for the n th coefficient of the polynomial: 

a n = (-lHA eiA% • ( A - 5 ) 



7 primitive 

Biy=n 



Appendix A. 2. Cancelations due to a combinatorial identity 

A pseudo orbit expansion for quantum graphs appeared in [M] . The formula for 
evaluating the pseudo orbits contributions was given there, but the coefficients were not 
computed explicitly We compute those coefficients and show that many cancelations 
occur due to a combinatorial identity, leading again to the formula (1A.5j) . 



In the following we will relate pseudo orbits to partitions of an integer. A 
partition of some positive integer, n, is an ordered set P = (mi, m.2, . . . , m n ) such 
that X^j=i 3 m j = n - m other words, it is a way to present n as a sum of all positive 
integers smaller than n (each can be taken once, more than once, or not at all). The 
set of all partitions of n is denoted by V n . 

We now relate a partition of n to a certain class of pseudo orbits. Let P = 
(m 1 ,m 2 ,...,m n ) G V n and denote by T P the set of all pseudo orbits, 7, such that 
for all 1 < j < n, 7 contains mj periodic orbits of topological length j. Hence for 
7 G Tp we have = n. 

We begin the proof by employing Newton's relations which connect coefficients of 
the characteristic polynomial a n to traces of powers of U (k) (see e.g. [20]). One form 
for these identities is 

1 n— 1 

a n = --J2^{U n ~ 3 (k))a j , (A.6) 

U 3=0 

with a = 1. We use this recursive relation to give a closed expression for the coefficients 
in terms of the traces. 

a = Y — (-IlY 1 J_ ( _TA m2 . . J_ (_TA mn (A 7) 
V 1/ m 2 ! V 2 / "' mj \ nj 

where Tj := tr (W (k)). A similar expression appears in [20l[3l]. In what follows we show 
that (TAJD can be expanded as a sum over pseudo orbits and that careful book-keeping 



allows a significant number of cancelations. 

The traces Tj can be presented as sums over periodic orbits on the graph as in 
(1A.2I) . Substituting in (1 A. 71) we obtain 



1 ( Tj 



^ n m 

PeVn 3=1 3 

En 



Pe^n 3=1 I J \ -rl s 7 =i 7 
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PeV n j=l { r \y\B^=j 7 / J 

One may expand the second sum and the product in the expression above. This results 
in a sum, each of whose terms corresponds to a product over periodic orbits. Such a 
product can be identified with a pseudo orbit. More specifically, 

a n = Cj A je ikh , ( A -9) 

; y\Bx / =n 

where is a coefficient that consists of combinatorial factors and repetition numbers. 
There are many cancelations in this sum and we will show that eventually only primitive 
pseudo orbits remain. 

Let us fix some primitive periodic orbit, n. Consider all pseudo orbits 7 with 
= n for some fixed n, which contain only copies of 7r and its possible repetitions. 
For example, such a pseudo orbit might be 7 = {5ir, 2tt 3 }, which means that 7 contains 
seven periodic orbits; five copies of tt and two copies of a periodic orbit where tt is 
repeated three times. For all such pseudo orbits = n is a multiple of B n , and we 
denote q := B^/B n . Summing the terms in ( 1A.9I) over this collection of pseudo orbits 
and reading the expression for from (1A.8j) . we obtain 



QdV q \r=l r 




where the sum is over partitions of q, Q = (ri\, . . . , n q ). To obtain (lA.lOj) observe that; 

(i) We only sum over pseudo orbits belonging to n-partitions in (I A . 8 [) with rrij = for 
all j values which are not multiple of B n . So we may enumerate these by considering 
the partitions of q = B^/B^ where n r = m r B^ ■ 

(ii) Each partition of q determines a unique pseudo orbit constructed only from tt. 
Denoting a g-partition by Q = (ni, . . . ,n q ), the repetition numbers of the orbits 
within our pseudo orbit are just the indices of the n r 's. 

(iii) All these pseudo orbits are weighted by the same factor, (A^e^y. 

Lemma [1] will show that the sum in (lA.lOj) is zero for all q > 1. This implies that 
the only surviving contribution comes from the pseudo orbit that contains 7r alone. 
Furthermore, we can extend this cancelation argument as follows. Fix n and let tt\ and 
7r 2 be two primitive periodic orbits, such that n = qiB ni + q^B^. We now examine all 
pseudo orbits with n bonds which contain only copies of 717, 712, and their repetitions. 



Collecting the contributions of all such orbits from flA.8j) . we arrive at 



Qi£V Ql \r=l r- \ / J Q 2&Vq2 \ T=1 r- \ / J 

Applying lemma [1] again, the sum above vanishes if either q\ > 1 or q 2 > 1 and hence 
the only contribution that does not vanish is for the primitive pseudo orbit {vr 1 ,7r 2 }. 
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Iterating this argument we conclude that all contributions from non-primitive pseudo 
orbits cancel. Hence 

a n = (-IH-V^ ( A - n ) 

7 primitive 



which agrees with (1A.5[) obtained using the first canceling mechanism. 



Lemma 1. For all q G N with q > 1, 



Proof. The expression we wish to vanish is equal to the coefficient multiplying x q in the 
following generating function, 

oo oo I / x r \ nr °° / X r 

n£^(-vj =ri ex p(-7 

r=l n r =0 




□ 



Appendix A. 3. Further diagrammatic cancelations 



Comparing the results of the two canceling mechanisms f ]A.5j) and (lA.llj) with theorem 



[U we see the difference is between a sum over all primitive pseudo orbits and a sum 
only over the irreducible ones. More specifically, in theorem [1] we do not sum over 
primitive pseudo orbits which contain the same directed bond more than once. It is 
possible to show that the contribution of all such periodic orbits cancel each other using 
a diagrammatic approach, as was shown very nicely in pp. One should note that in [1] 
the diagrammatic cancelations are presented assuming each vertex scattering matrix, 
o~( v \ has all its off diagonal entries equal. However, their diagrammatic cancelations 
which can be used to transform ( lA.llj) into ( l30i) do not require this assumption. We 
also note that a similar cancelation mechanism appears in [8]. However, the derivation 
there is done for the dynamical zeta function and the cancelations occur between orbits 
which pass through the same cell in the Poincare surface of section. 
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